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O : 1 INTRODUCTION 

Although the presently most promising approach to the direct detection of gravita- 
tional radiation would seem to be provided by the use of appropriate electromagnetic 
field configurations (e.g. that of and optical interferometer) the original, and until 
^ . now most highly developed method has been based on the use of an elastic solid con- 
figuration, traditionally a cylindrical bar of the kind originally introduced by Weber. 
A highly simplified but essentially adequate description of the interaction of gravita- 
tional radiation with such an apparatus was given by Weber himself in his seminal 
1960 paper [1]. For a more detailed analysis, and for the treatment of more general 
configurations such as that of the earth as a whole, a system of wave equations gov- 
erning the interaction of weak gravitational radiation with an elastic solid was derived 
soon after by Dyson [2], Papapetrou [3], and others, subject to the limitation that 
non-linearities due to self gravitation of the solid medium are neglected, as is entirely 
justifiable in a terrestrial context. 

While the final outcome of the present article will be a rederivation of the weak 
field wave equations to which we have just referred, we shall nevertheless use a quite 
different approach to that of Dyson and Papapetrou. Instead of working throughout 
in an only approximately self consistent linearised scheme, as they did, we shall first 
set up the exactly self consistent fully non-linear theory of the interaction of a gravi- 
tational field with an elastic solid in accordance with Einsein's theory. The general. 
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non- linear theory is in any case needed for application to the more exotic context of 
neutron star deformations, as was discussed e.g. by Carter and Quintana [4]. Start- 
ing from this mathematically sound basis, we shall then proceed to the derivation of 
the weak field limit in two successive stages of approximation: we shall first impose 
the restriction that the gravitational radiation be weak, even though the unperturbed 
background field may still be strong (as in the case of a neutron star); finally we shall 
impose the condition that the background field also be weak (as in the case of the 
earth) so as to obtain the Dyson Papapetrou equations. 

Unlike the treatments used by other authors such as those of the schools of Eringen 
(see e.g. Maughin [5]) or of Souriau (who was the first to set up the fully non-linear 
elasticity theory used here [6]) our present treatment will be based on the use of 
convected differentials, a powerful technique for the analysis of material media (elastic 
or otherwise) that generalises the more restricted convected differentiation procedure 
first introduced by Oldroyd [7]. In providing a self-contained introduction to the 
general concept of convected differentals and differentiation, and to their application 
to the theory of an elastic solid and its interaction with gravitation, the present article 
condenses several previous publications describing work carried out in collaboration 
with H. Quintana [8], [9], [10], [11], [12]. This article thus serves to update an 
earlier introductory survey by Ehlers [13] . 

2 KINEMATICS OF A MATERIAL MEDIUM: MATERIAL REPRESENTATION 

In Newtonian theory a material medium is usually visualised as a three dimen- 
sional manifold whose configuration at a given instant is specified by a non-singular 
mapping into three dimensional Euclidean space, the motion of the medium being 
given by the time variation of the maping. However in General Relativity theory, 
where there is in general no canonically preferred set of three-dimensional sections 
(Euclidean or otherwise) of the fundamental four-dimensional spacetime manifold A4 
say, it is more convenient and natural to proceed the other way about. The medium 
itself is still to be conceived in the abstract as a three-dimensional manifold, X say, 
(whose points represent idealised particles of the material) but its motion can be 
specified by a necessarily degenerate (four to three dimensional) mapping, V say, of 
the world-tube traversed by the matter in M. onto the manifold X, the world line of 
each idealised particle of the medium being projected under V onto the corresponding 
point in X. 

We shall use Greek indices /x, i/, ... to specify ordinary spacetime tensors onAi, and 
capital Roman indices A, B, ... to specify tensors on the manifold X representing the 
medium. The spacetime indices i^, ... may be thought of in the traditional manner 
as specifying tensor components defined in terms of a local system of coordinates x'^ 
{p = 0, 1, 2, 3) on Ai while similarly the material indices A, B, ... may be thought of as 
defined in terms of local coordinates (a=1,2,3) on X. However, following Penrose 
[14], it wil also be comvenient to interpret the index symbols in the abstract manner 
as merely an indication of the tensorial (or more general) character of the quantities 
concerned (rank, co/contravariant quality, and in the present case association with A4 
or X) rather than as integers specifying concrete components. Thus with the indices 
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interpreted in this loose sense, the statement that a point with coordinates in 
X is the image under P of a point with coordinates in M. may legitimately be 
expressed in a simple and natural manner by 

= V{x^'} (2.1) 

whereas with a strict traditional interpretation of A and such an equation would be 
nonsense. 

The operation V of projection from Ai onto X evidently induces a corresponding 
projection, which we shall also denote by V, from tangent vectors at any given point x^^ 
in M. onto tangent vectors at the corresponding point 'P{x^} in X. In an obviously 
natural notation we shall denote the projected image of a tangent vector by the 
same basic symbol, distinguishing the image from the original vector only by the 
appropriate change from Roman to Greek of the index symbol, i.e. for any spacetime 
tangent vector we shall set 

e = v{e} ■ (2.2) 

Due to the degeneracy of V there is in general no corresponding induced projection 
operation for covectors except in the particular case of a covector ck^ that is orthogonal 
to the congruence of worldlines in the sense that 

uf'a^^O, (2.3) 

where the vector is tangent to the worldline at the spacetime point in question. 
In this particular case the projection will induce a corresponding covector on X, for 
which as before we shall use the same basic symbol, i.e. we shall write 

«A = V{a^} , (2.4) 

and in this case the operation is reversible, i.e. we shall have a bijection 

CK/u ^ ttA (2.5) 

between orthogonal covectors at in M. and covectors at V{x^} in X. We can use 
the (pseudo-) metric tensor g^^, on the spacetime manifold A4 (with signature — 
and units such that c = 1) to define a corresponding bijection for tangent vectors by 
restricting our attention to tangent vectors on M. that are orthogonal to the world 
lines in the metric sense, i.e. 

Uf,e = 0, u^=g^.u''. (2.6) 

Subject to this restriction the projection (2.1) will have the same reversibility property 
as (2.4), i.e. 

e ^ e (2.7) 

in the same sense as (2.5). 
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The natural 1-1 correspondence that has just been defined between vectors or 
covectors orthogonal to the worldlines in spacetime and vectors or covectors in the 
material medium can evidently be extended directly to general orthogonal tensors as 
defined in terms of tensor products of orthogonal vectors and covectors. This corre- 
spondence is of vital importance for setting up any physical theory of the mechanical 
behaviour of the medium, since its spacetime evolution must be described in terms of 
geometrical quantities (tensors etc.) defined in spacetime, whereas its intrinsic prop- 
erties can only be described in terms of quantities described directly in terms of the 
manifold X representing the material medium (a requirement commonly dignified as 
the "principle of material objectivity"). 

Since there is no guarantee that all physically relevant spacetime tensors will turn 
out to be automatically orthogonal to the worldlines, it is important to remark that 
they can always be canonically decomposed into a set of orthogonal tensors of equal or 
lower order by using the orthogonal projection tensor 7^. In terms of the normalised 
unit tangent vector u^^ and the corresponding covector defined by 

uf'u^ = -1 (2.8) 

this orthogonal projection operator is defined in terms of the unit tensor gf^ by 

7^^ = < + u^'u, . (2.9) 

In the particular case of an ordinary tangent vector the natural decomposition into 
a set consisting of an orthogonal part ±v^ and a scalar is given by 

^^v^" + vK'^ , (2.10) 

where 

^v^' = ^^^v'', v^^ = -u^v^' , (2.11) 

while for a general covector a;^ we similarly have 

i^i^=±i^fi + i^\\u^, (2.12) 

where 

-L^M = T^'^i^ ) '^11 = -u'^'^fi ■ (2.13) 

This allows us to represent them in terms of sets of corresponding geometrical quan- 
tities defined at the corresponding point in the material medium in the form 

yf^ ^ {^v^,v^ (2.14), 

and 

^ {±^A,(^f,}, (2.15) 

the extension to general tensors of higher order being straightforeward albeit cum- 
bersome. 
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3 KINEMATICS OF A MATERIAL MEDIUM: CONVECTED DIFFERENTIALS 



The concept of a material representation of spacetime tensors in terms of sets 
of tensors defined on the medium space (as defined by (2.14) and (2.15) ) gives 
rise naturally to the concept of material variation defined as the difference between 
the material representations in any Lagrangian (i.e. worldline preserving) variation 
between different configurations, which may either arise from a mapping between dif- 
ferent physically conceivable evolutions of the medium or else merely from a time 
displacement in a single evolution. The spacetime tensor corresponding to the in- 
finitesimal material variation between nearby states of evolution of the medium will 
be referred to as the convected differential [12]. We shall use the symbol A to denote 
an ordinary Lagrangian differential and we shall use the notation d[ ] for a correspond- 
ing convected differential. Thus in the case of a vector and a covector respectively, 
(2.14) and (2.15) give rise to the correspondences 

d[v^] ^ {Ajv'',Av^, (3.1) 

and 

d[uj] ^ {A_La;A, Aa;||}. (3.2) 

To evaluate these expressions we use the facrthat the world line preserving variation 
can only change the magnitude but not the direction of the unit tangent vector u^^: 
explicitly 

Au" = ^u^u^u^Ap^ (3.3) 

where we use the abbreviation 

A,^ = Agp^ (3.4) 
for the Lagrangian variation of the metric. In the case of a vector we obtain 

Avil = -u^.Av" - v^Au^ (3.5) 



A^^;'^ = 7^A^;^ + v'' {u^'Au^ + u^Au") . (3.6) 

Since the bracketed quantity in the last expresion is automatically parallel to u'^ it 
does not contribute to the material projection, so we obtain 

A^v"- = ViAjv^"} = Vi-i'^v") (3.7) 

and hence 

dK] =7^At;'^ + w'^At;II 

= Av^" - u^v'^Au^ . (3.8) 

In the case of a covector we obtain 

Acc;|| = -vy Aujy - ujyAvy (3.9) 
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^Apj^i = Y^^oj^ + iOv{vyi^u^ + u^,/\vr) , (3.10) 

both terms in the last expression being automatically orthogonal to u^^. Thus using 

AauJa = r{Ajuj^} (3.11) 

we immediately obtain 

d[ujij] = A_La;^ + it^Aa;|| 

= Ao;^ + uji,u'^Aun . (3.12) 

The extension to a general tensor is now automatic: it suffices to add an appropriately 
analogous term for each extra index. Thus for a general mixed tensor T^ ;; one obtains 
the relation [12] between convected and Lagrangian differentials in the form 

dlTt-;;] = AT^tz + T%-- u^Au^ + ... 



-TP- ui'Aup - ... (3.13) 

where the Lagrangian differential of the covector is obtainable by substituting the 
formula (3.3) in the expression 

Au^ = g^uAu" + u'^A^^ . (3.14) 

An important particular case [8] covered by the general formula (3.13) is that in 
which instead of making a comparison between nearby but different states of material 
motion (as one needs to do in perturbation theory) one wishes to study time variations 
in a single given state of motion. This corresponds to the case in which the Lagrangian 
variation is simply given by Lie differentiation with respect to a time displacement 
vector field, C,^^ say, which we shall denote by C£, i-e. we set 

A = C£, (3.15) 

where the vector C,^ is an arbitrarily normalised tangent to the flow, which can there- 
fore be expressed in the form 

C^' = u^'AT, (3.16) 

where At is an arbitrary scalar field interpretable as representing the local value of the 
corresponding infinitesimal proper time displacement along the world lines. Although 
the various terms in (3.13) will involve derivatives of C,^ , the intrinsic nature of the 
material variation will ensure that the convected differential will only depend on the 
scalar value of the displacement dr at the point under consideration, so that it will 
take the form 

d[T^.::] = [T^.::]dT (3.17) 
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where the tensor [T^;; ]' so defined is what we refer to as the convected derivative . 
By direct substitution of (3.15) and (3.16) in (3.13) one can check that the terms 
involving gradients of dr do indeed cancel out, so that one recovers the original formula 
[8] for the convected derivative, namely 

[n.-.-.r = n.;: + t%:_ + v.)ii^ + ... 

-TP;::XuP + Vp)u''-... (3.18) 

where is the usual metric covariant differentiation operator and where we use a 
simple dot withoutsquare brackets to denote covariant differentiation with respect to 
the propert time, i.e. ' = u'^'Vp. In the particular case of the unit tangent vector 
itself, the dot operation gives the acceleration vector, whose covariant form is also 
expressible as the Lie derivative: 

= u£ui^ . (3.19) 

Just as (3.13) is a generalisation of the earlier formula (3.18), so also (3.18) is 
itself a generalisation of a previuos formula given by Oldroyd [7] for the particular 
case of tensors entirely orthogonal to the world lines, for which, as pointed out by 
Ehlers [13], the convected derivative reduces to the orthogonal projection of the Lie 
derivative. An important example is the divergence tensor O^iy of the material flow, 
as defined by the decomposition 

V^Ujy = Op,i, + UJ^y - Up,Uy , (3.20) 

where the vorticity tensor cupi, is antisymmetric and 9^1, is symmetric. It is related to 
the strain tensor j^^, (i.e. the covariant version of the orthogonal projection tensor 
(2.19) ) by 

Op.^^M. (3.21) 

4 MECHANICS OF A PERFECT ELASTIC MEDIUM 

The convected differential and derivative that were described in the previous sec- 
tion are potentially useful for the kinetic analysis of any kind of material medium. One 
of the simplest applications is to the theory of a medium whos behaviour satisfies the 
following (by now standard) criterion of perfect elasticity. A perfect elastic medium 
can be characterised succinctly by the condition that its energy - momentum tensor is 
a material function of the metric tensor with respect to the flow fleld specifled by its 
timelike eigenvector. This means that the energy momentum tensor takes the form 

T'^^ = pui^u" + P"'' (4.1) 

where the pressure tensor satisfies the orthogonality condition 

pf^'^m, = , (4.2) 
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so that the material representation of T^^ is expressible as 

r/^- ^ {P^^,0,p}, (4.3) 
and this representation must be a function of the corresponding representation 

9f.u ^ {tab, 0,-1}. (4.4) 

Thus at each point in the three dimensional manifold X representing the medium 
there are well defined functions determining the pressure components P^^ and also 
the density p in terms of the components 7ab- There is however a restriction that 
prevents these functions from all being chosen arbitrarily, namely the local energy - 
momentum conservation law 

V^T^" = 0, (4.5) 

which has four independent components, whereas the acceleration of the flow has only 
three independent degrees of freedom. In order to avoid having an overdetermined 
system of equations of motion, one must require thqt the component of (4.5) along the 
flow (i.e. the conservation of rest frame energy as distinct from momentum) should 
be satisfled as an identity. The remaining independent equations are given by 

T^VpT^^^O, (4.6) 

which is equivalent to the equations of motion 

pu^' = -^i^^^pp^ . (4.7) 

The equation that must be satisfled identically is 

ti, V^T'^" = , (4.8) 

which is equivalent to 

p=_(p^M- + pM-)^^^. (4.9) 

Now it follows from (3.21) that this last can be expresssed in terms of convected 
derivatives as 

\p] = -\^p^^^^P^^)Yi^,]^ (4.10) 
which means that the convected variations allong the world lines must satisfy 

d\p\ =-\{pi^''^P^'')dYi^..\. (4.11) 

and hence that the corresponding variation of the material projections in X must 
satisfy 

^^P = -Mp7^''+^^^'')^7ab. (4.12) 
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This will be satisfied automatically for an arbitrary equation of state p = p{7ab} if 
and only if the corresponding equations for the six algebraically independant pressure 
components are specified by 



By carrying out a second partial difi'erentiation of the single equation of state 
function for p with respect to the strain 7ab we deduce that the material variation of 
the pressure tensor will be given by 

^pAB ^ _1 ^^ABCD ^ pAB^CD^^^^^ ^ (4_;L4) 

where the elasticity tensor, whose material projection is defined by 

^ABCD ^ _ pAB CD 

will obey the symmetry conditions 

^Hvpa _ ^{nv){pa) _ ^panv (4-16) 

as well as the orthogonality requirement 

E^'^P'^Ua^O. (4.17) 

A familiar special case of a perfectly elastic medium is that of an ordinary perfect 
fluid, which can be defined in the present context by the condition that its density be 
a function only of the determinant I7I of the material projection of the strain tensor. 



I.e. 



p = p{|7|}. (4.18) 



It follows from (4.13) that the pressure tensor will then take the purely isotropic 
form 

pM^ = P7^^ , (4.19) 



with the pressure scalar given by 



P = -2|7l^, (4.20) 



while the elasticity tensor will be given in terms of the bulk modulus 

/5 = -2|7l^. (4.21) 
9 



by the formula 



(4.22) 



5 SMALL GRAVITATIONAL PERTURBATIONS OF AN ELASTIC MEDIUM 

Having seen how to set up a system of exactly self - consistent but non - linear 
equations governing a perfect elastic medium in the framework of General Relativity, 
we are now ready to derive the linearised wave equation governing small perturba- 
tions relative to a known background, such as might be induced by weak incoming 
gravitational radiation. 

It is usually convenient to think of the perturbations as being determined in terms 
of a vector field that specifies the infinitesimal displacement of the worldlines rela- 
tive to their positions in the known background space. Of course such a displacement 
is entirely gauge dependent and can always be reduced to zero by the use of an ap- 
propriate Lagrangian (worldline dragging) mapping of the perturbed space onto the 
background, butit is often convenient to fix the gauge by a more purely geometric 
requirement such as the preservation of a harmonic coordinate system, which is the 
generalisation of the usual newtonian procedure of defining the displacements relative 
to the fixed Euclidean structure of space. We shall use the symbol S to denote the 
Eulerian variation of any quantity as specified by any such geometric prescription 
for the mapping of the perturbed spacetime onto the unperturbed background. The 
difference between the Lagrangian variations denoted by A and the Eulerian varia- 
tions denoted by S is given by Lie differentiation with respect to the corresponding 
infinitesimal displacement field i.e. 

A = S + ^£. (5.1) 

Thus in particular the Lagrangian variation of the metric tensor is given by 

A^^ = V + 2V(^e^) , (5.2) 

where we use the usual notation 

hij,iy = 5g^y (5.3) 

for the Eulerian variation of the metric arising from the gravitational waves under 
consideration. 

As far as quantities characterising the material medium are concerned, it is easier 
to work with Lagrangian than Eulerian variations, since the former are related directly 
(via (3.13)) to convected variations and hence to the material variations that are 
governed directly by the equations of state. In the case of orthogonal covariant tensors, 
including the special case of scalars, we see from (3.13) that the Lagrangian variation 
is given directly by the convected variation, and hence from (4.12) we find that the 
Lagrangian variation of the density is given in terms of that of the metric by 

Ap=-ip|/'"^Ap,, (5.4) 
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where for future convenience we introduce the abbreviation 



yP-^ = + p-^pP'' . (5.5) 

For a contravariant tensor however, (3.13) introduces extra terms, so that for the 
dependence of the Eulerian variation of the pressure tensor in terms of the metric 
(4.14) gives rise to the formula 

^pM^ = _i(^M^P'- + py-^Y" - 4PP^^u''^u'')Ap^ . (5.6) 

Using these results together with the expression (3.3) for the lagrangian variation 
of the flow field itself, we are now in a position to work out the perturbed equations 
of motion by taking the variation of the exact equations of motion (4.7). It is evidently 
most convenient to start from the Lagrangian variation, i.e. 

A(pu^ + 7^;VpP^^)=0, (5.7) 

which works out explicitly as 



= [P^'^'if - \ui'P'"' - 2yl'^(; ^) {QP^ + ujP^)u'' + py^pUpyyu''^ Ka (5.8) 
using the abbreviation 

for the relative strain tensor, and 

for the modified elasticity tensor first introduced in a classical context by Hadamard. 
When the Lagrangian metric perturbation A^,^ is evaluated by use of (5.2), and 
the corresponding perturbation of the affine connection components T^^^^ is evaluated 
using the corresponding substitution 

AL'^,, = V(,A'^,) - ^V^A,, 

= V(,V,)e^ + V^^W,^ - \V^K, - eR%,)p (5.11) 
where the Riemann tensor is defined by 

(V^V, - V,V^)ep = Rf^upaC , (5.12) 
11 



then we see that the basic perturbation equation (5.8) takes the form of a hyperbohv 
wave equation for the displacement vector when is given. The characteristic 
cones and the corresponding sound speeds can be worked out directly from (4.7) 
by considering discontinuities [10] without any need of the full set of perturbation 
equations given here. 

In order to have a complete system of equations governing the interaction of weak 
gravitational radiation with an elastic medium we need an additional wave equation 
governing the gravitational perturbation h^,^. This is obtainable by taking the ap- 
propriate perturbations of the Einstein gravitational equations, which are expressible 
(in units with G=l) by 

Rf" = SttT^"" (5.13) 
where the Ricci tensor is defined by 

Rf^t^ = R^pf (5.14) 
and where we use the notation "for the partial trace subtraction operation defined by 

R^'^ = R^'- - Igi^- Rp . (5.15) 

In our work up to this stage we have been concentrating on material aspects, so 
that Lagrangian variations have given the simplest formulae. However the advantages 
of being able to use more general Eulerian variations become apparent now that we 
come to the properly gravitational aspect, since it is well known that the perturbed 
Einstein equations in the Eulerian form 

- SttT'^^) = (5.16) 

can be greatly simplified by the imposition of the harmonic gauge condition 

V^/i'^^ = 0. (5.17) 

Under these conditions most of the terms drop out and one is left with an equation 
of the form 

□ /i/^^ = -167r(5r''^ (5.18) 
where the relevant wave operator is defined by 

-\Rf>{h'^--lg^^hP) (5.19) 

using the standard notation 

^'^p. = R'^'p. - 2^1"pK^] - k9l^R\) (5.20) 
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(with square brackets denoting antisymmetrisation) for Weyl's trace free conformal 
tensor. The Eulerian variation of the energy-momentum tensor is obtainable using 
(5.1) in terms of a Lie derivative and a more easily evaluable Lagrangian variation. 
Thus starting from (4.1) and using (3.1), (5.4), and (5.6) we obtain finally 



+2TP^i'VpC^ - C^pTf"" (5.21) 

where, following Friedman and Schutz [14] we have constructed a generalised (non- 
orthogonal) elasticity tensor with the same symmetry properties as those of the ordi- 
nary elasticity tensor (4.16) according to the prescription 

u'^u'^ . (5.22) 

The coupled system of equations (5.8) and (5.9) simplifies considerably when 
not only the perturbations but also the background gravitational field is weak, as is 
the case in terrestrial (as opposed to neutron star) applications. In such cases we may 
suppose that there is a small dimensionless parameter, e, loosely interpretable as an 
upper bound not only on the order of magnitude of the gravitatinal wave perturbations 
hp,iyi e^jy etc., but also of the deviations of the metric from the flat Minkowski form. 
It then follows from the Einstein equation that the density p must also be of linear 
order in e as the latter tends to zero, while (by the virial theorem) the pressure in a 
self gravitating system is of even higher order, tending to zero even when divided by 
e, i.e. in standard notation 

P''" = o{e} . (5.23) 

(In our original version [11] a printer's error substituted in place of o throughout.) 
Since the unperturbed energy-momentum tensor will be at most of linear order in e 
its perturbation will be of higher order, i.e. 

ST"'' = o{e} (5.24) 

so that the gravitational wave equation (5.18) will to lowest order be of simple 
Dalembertian form , i.e. 

WpW'h'"' = o{e} . (5.25) 

It follows that in addition to the harmonic gauge condition we can to this order make 
the further simplifications 

= o{e} , /i"^ = o{e} (5.26) 

and to the same order the wave equation for the displacement reduces to the form 
w^u'^V.V^e - p-'V^iE^'-'^Vp^a + ^ V)) = o{e} (5.27) 
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in agreement with calculations of Dyson [2] and Papapetrou [3] . 

As a simple practical application Dyson considered the case of an elastic medium 
that is isotropic, as it will be in the case of a typical metal when considered on 
scales large compared with the microscopic crystalline domains. In such a case (in 
consequence of (5.23)) the elasticity tensor will be of the form 

Et^i^pa = 13^1^'^ + 2|I{-f^'^P-f''> - ^^^''Y") (5.28) 

where (5 is the bulk modulus and /i is the rigidity modulus. Now under these conditions 
it follows from (5.26) that the only gravitational term in (5.28) reduces to the form 

V. {E^'^'P'^hp,) = 2/^^^^ V./x , (5.29) 

which shows, as pointed out by Dyson, that the gravitational waves couple with the 
elastic displacement only via non-uniformities of the rigidity. (In the case of a tradi- 
tional Weber bar detector the relevant non-uniformity is provided by the discontinuity 
at the surface of the cylinder.) 

REFERENCES 

[I] Weber, J., Phys. Rev. 117 (1960) 306. 
[2] Dyson, F.J., Astroph. J. 156 (1969) 529. 

[3] Papapetrou, A. Ann. Inst. H. P oincare A16 (1972) 63. 

[4] Carter B., and Quintana, H., Astroph. J. 202 (1975) 54. 

[5] Maugin, G.A., J. Math. Phys. 19 (1978) 1198. 

[6] Souriau, J.M., Geometric et Relativite (Herman, Paris, 1965). 

[7] Oldroyd, J.G., Proc. Roy. Soc. A270 (1970) 103. 

[8] Carter, B., and Quintana, H., Proc. Roy. Soc. A331, (1972) 57. 

[9] Carter, B., Commun. Math. Phys. 30 (1973) 261. 

[10] Carter, B., Phys. Rev. D7 (1973) 1590. 

[II] Carter, B., and Quintana, H., Phys. Rev. D16 (1977) 2928. 
[12] Carter, B., Proc. Roy. Soc. A372 (1980) 169. 

[13] Ehlers, J., in Israel, W., (ed) Relativity, Astrophysics and Cosmology (Reidel, 
Dordrecht, 1970) 89. 

[14] Penrose, R., in DeWitt, C, and Wheeler, J. A. (eds) Battelle Rencontre (Gordon 
and Breach, New York, 1968) 

[15] Friedman, J.L., and Schutz, B.L., Astrophys. J. 200 (1975) 204. 



14 



